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Abstract 

The dynamics of a geometric measure of the quantum discord (GMQD) under decoherence is 
investigated. We show that the GMQD of a two-qubit state can be alternatively obtained through 
the singular values of a 3 x 4 matrix whose elements are the expectation values of Pauli matrices 
of the two qubits. By using Heisenberg picture, the analytic results of the GMQD is obtained 
for three typical kinds of the quantum decoherence channels. We compare the dynamics of the 
GMQD with that of the quantum discord and of entanglement. We show that a sudden change in 
the decay rate of the GMQD does not always imply that of the quantum discord, and vice versa. 
We also give a general analysis on the sudden change in behavior and find that at least for the Bell 
diagonal states, the sudden changes in decay rates of the GMQD and that of the quantum discord 
occur simultaneously. 
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I. INTRODUCTION 



Correlations of bipartite states, including classical and quantum parts, are of great im- 



portance and interest in quantum in 
quantify the quantum correlations [l. 



ormation theory. Quantum discord was proposed to 
2|. It was suggested that the quantum discord, rather 



than entanglement, is responsible for the efficiency of a quantum computer, which is con- 
firmed both theoretically [3] and experimentally j4|. A great deal of efforts has been devoted 
into the study of quantum discord [l-27|. Despite this, it is not easy to obtain analytical 
results of the quantum discord since the optimization procedure involved is unreachable for 
arbitrary bipartite states up to now. Even for two-qubit systems, the analytic results are 
only known for a few cases jsUl, and a general method still lacks. To avoid this difficulty 
and obtain an analytic analysis, alternative approaches are needed, among which is the ge- 
ometric measure of quantum discord (GMQD). Despite not reflected in the present work, 
another advantage of the GMQD is that it could potentially supply a way to put various 
correlations on an equal footing since the geometric measure of other kinds of correlations 
can be deflned in the same manner but with different sets of zero-correlation states. It 
is remarkable that a unifled view of correlations has been established through the relative 
entropy measure of corre 



sure of the entanglement 



ations in Ref. 



28 



20|. The GMQD, similar to the geometry mea- 



29|, is deflned as the nearest distance between the given state 



and the set of zero-discord states. In the present work, we use the Hilbert- Schmidt norm 
as the distance of two quantum states, because for two-qubit systems the minimization of 
the Hilbert-Schmidt distance over the set of zero-discord states was resolved analytically in 



Ref. 



19| . However, the quantum discord is based on the Von Neumann entropy while the 



GMQD is based on the geometric distance, their behaviors may be different. This motivates 
us to consider the behaviors of the GMQD under decoherence, and compare it with the 
quantum discord. 

Due to the inevitable interaction with environment, the dynamics of the quantum discord 



under decoherence is of great importance. It has received some investigations 12 



17j . and 

lat 



15| 



was experimentally investigated in an all-optical setup most recently [18]. It was found t 
the quantum discord may decay in an asymptotic way under Markovian environment 
and vanish only at some time points under non-Markovian environment JigI . [l^. It can be 
understood by the facts that the subset of the zero-discord states has measure zero and is 



2 



nowhere dense 



22|. While the entanglement suffers form sudden death |30M32| . because the 



set of separable states occupies finite volume 



33l435ll. Besides, in some situations, the decay 

nn 



rates of the quantum discord may be discontinuous |12|. Il4|. This is a novel phenomena and 



was observed in the recent experiment 



18| . Notice for a tripartite system ABC in pure 



states, the quantum discord of AB and the entanglement of formation (EoF) of AC are 



connected through a monogamy relation [36|, l37| , so studying the dynamics of the quantum 



discord will also be helpful to the understanding of the dynamics of EoF. 

In the present work, we investigate the GMQD under decoherence channels and get the 
analytical results. Under three typical quantum decoherence channels, we will show that 
the GMQD is monotonically non- decreasing with respect to the quantum discord. Yet, the 
quantum discord may keep constant while the GMQD decreases. We will show that in some 
cases the decay rates of the GMQD and of the quantum discord may suddenly change at the 
same time. However, a sudden change in the decay rate of the GMQD does not always imply 
a sudden change in the decay rate of the quantum discord, and vice versa. We demonstrate 
each case by instances. We also give a general analysis for the sudden change in the decay 
rate of the GMQD and that of the quantum discord, and show that at least for the Bell 
diagonal states, the sudden changes in decay rates of the GMQD and that of the quantum 
discord occur simultaneously. 

This paper is organized as follows. In Sec. [Tll we give a brief introduction of the quantum 
discord and the GMQD. Then we show that the GMQD of a two-qubit state is related to the 
singular values of a peculiar 3x4 matrix. In Sec. IIIIl we give a general method to obtain the 
GMQD under quantum decoherence channels, and get the analytic results for three typical 
kinds of quantum decoherence channel. We investigate the sudden change in the decay rate 
of the GMQD, and compare it with the case of the quantum discord. And we also give a 
general analysis on the disagreement of sudden change in decay rates between the GMQD 
and the quantum discord. Section [IV] is the conclusion and discussion. 



II. GEOMETRIC MEASURE OF QUANTUM DISCORD 



Given a quantum state p in a composite Hilbert space = T-La ®'Hb, the total amount 
of correlation is quantified by quantum mutual information 



Z{p) = H{pa) + H{pb)-H{p), 



(1) 



where H{p) = — Tr [plog2 p] is the von Neumann entropy and Pa{b) = '~^t^b(A)P is the reduced 
density matrix by tracing out system B[A). If we take the system A as the apparatus, the 
quantum discord is defined as follows [l, 2| 



VAip)=Iip)-CA{p), 



(2) 



which is the difference of the total amount of correlation X(p) and the classical correlation 
Ca{p)- Here the classical correlation is defined by 



Ca{p) = maxX(p|{Efc}), 



(3) 



where X{p\{Ek}) is a variant of quantum mutual information based on a given measurement 
basis {Ek} on system A as follows 



I{p\{Ek}) = H{pb) - Y,PkH{pB\k)- 



(4) 



PB\k = Tr^[(_E'j!c ® t.)p\/pk is the postmeasurement state of B after obtaining outcome k on 
A with the probability pk = Tr[(ii^j!c ® l)p]- {Ek} is a set of one- dimensional projectors on 
I-Lai and 1 is the 2x2 identity operator. 

In Ref. 19|, Dakic et al. proposed a geometric measure of quantum discord defined by 



mm p-x 

xef^o 



(5) 



where VLq denotes the set of zero-discord states and | |X| p := Ti[X'^X) is the Hilbert-Schmidt 
norm. The subscript A of D\ implies that the measurement is taken on the system A. For 
two-qubit systems, a zero-discord state is of the form x = ® pi + P2\i^2) {"^21 ® P2 

with lipi) and \1IJ2) two arbitrary orthogonal states. And a general state can be written in 



Bloch representation 
P = 



3 3 

i i,j=l 



(6) 



with Xi, yi, and Rij real parameters, and o"j=i,2,3 Pauli matrices. Then an explicit expression 
of the GMQD is obtained as (loj: 



1 



^^(p) = T(lkll'+ll^ir-^max) 



(7) 



where x = {xi,X2, x^)'^, R is the matrix with elements Rij, and fcmax is the largest eigenvalue 
of matrix K = xx^ + RR^. 
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Now, we introduce an alternative form which will be convenient when we consider the 
evolution of the GMQD under decoherence. First, we introduce a matric TZ defined by 

(8) 

and another 3x4 matric TZ' obtained through deleting the first row of R, i.e., TZ' = (x, R). 
Here TZ is just the expectation matric with the elements TZij = Tr[pcrj(8)crj] for i,j = 0, 1, 2, 3, 
and (Jo = 1 is defined. The definition of TZ' leads to K = TZ'{TZ'y . After singular value 
decomposition, we have TZ' = UAV'^, where U and V are 3x3 and 4x4 orthogonal matrices, 
and A has only diagonal elements Aij = Xi6ij with A, the so-called singular values of the 
matrix TZ'. Then the eigenvalues of the matrix K can be expressed as A^. Considering 
+ ||-R|P = Tri^, we get an alternative compact form of Z)^(p): 

(9) 

where the summation and maximization are taken over all the non-zero singular values A^ of 
TZ'. This alternative form will be convenient when we consider the evolution of the GMQD 
under decoherence. 

III. GEOMETRIC MEASURE OF QUANTUM DISCORD UNDER QUANTUM 
DECOHERENCE CHANNELS 

A quantum channel can be described in the Kraus representation 

£ip) = J2K,pKl (10) 

where are Kraus operators satisfying "^^Kj^K^ = 1. As we discussed in the previous 
section, to obtain the GMQD, we need to know the expectation values of the Pauli matrices 
of the two qubits for the state S{p). So we turn to the Heisenberg picture to describe 
quantum channels via the map [40|] 

£^iA) = J2K^K, (11) 

with A an arbitrary observable. Then the expectation value of A can be obtained through 
{A) = Tr [y4£'(p)] = Tr [£'''(yl)p]. Because an arbitrary Hermitian operator on can be 
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TZ 



1 

X R 



DM 



1 
4 




max A? 

k " 



expressed hj A = ^^=o^«^« with G M, then a quantum channel for a qubit can be 
characterized by the transmission matrix M defined through 



(12) 



Since Tt[£'^ {cXi) p] = MjjTr[crjp], Mjj actually describes the transformation of the polar- 
ized vector Pi = Tr[(7jp]. 

Now we consider the case of two qubits under local decoherence channels, i.e., p = 
[£a ® Sb]{po)- To obtain the GMQD of the output state p through the channel, we need to 
get the expectation matrix TZ. With the Heisenberg picture, we have 



(13) 



where TZo is the expectation matrix under po, i-e., {'R-o)ij = Tr(o"j ® cTjPo), and Ma{b) is 
the transformation matrix characterizing the quantum channel Sa(b)- So we obtain TZ = 

For simplicity, we assume E"^ and £^ be identical, hereafter. Next, we consider three 
typical kinds of decoherence channels: the amplitude damping channel (ADC), the phase 



damping channel (PDC), and the depo 



set of Kraus operators respectively 



41 



arizing channel (DPC). They are described by the 



42|: 



{yi|o)(o| + |i)(i|, v^|l)(0|} 

/^l, v^|0)(0|, v^|l)(l|}. 



^ADC 
^PDC 



(14) 
(15) 
(16) 



with s = 



setup [41 



— p. Here the real parameter p G [0, 1] may be time-dependent in some realistic 



42| . For instance, for the PDC, the parameter s may be like exp(— 7t) with 7 the 
rate of damping. 

From Eqs. (1121) . (IT^ . ( JTSl) . and (fT6l) . the transmission matrix M of each channel can be 
got through the transformation of the Pauli matrices in the Heisenberg picture |4o| as 



(17) 
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For simplicity, here we first take as the input states of two-qubit system the Bell diagonal 



states 15 



12| 



P 



(18) 



which includes the Werner states (|ci| = |c2| = \c^\ = c) and Bell states (|ci| = |c2| = \cz\ 
1). This state is physical if the vector (01,02,03) belongs to the tetrahedron defined by the 
set of the vertices (—1, —1, —1), (—1, 1, 1) 



be described by the following conditions 



1, — 1, 1) and (1, 1, —1) 43|. This restriction can 



43|: 



G [-3,1], 



i=l 



Ci - Cj - Ck e [-3, 1] for j ^k. 



(19) 



For states (fTSj) . IZq = diag{l, Ci, C2, C3} is of diagonal form. From the relation 71 
M7^oM^, we get 7^ under the ADC, the PDC, the DPC respectively: 



"7^ ADC 
'T^DPC 



10 -p 
CiS 

C2S 
—p c^s^ + p^ 
diag {l,CiS^C2S^C3} , 



diag {l,cis^C2S^C3S^} . 



(20) 

(21) 
(22) 



TZ' is obtained by deleting the first row of the matrix TZ, for ADC, PDC, DPC respec- 
tively. Calculating the singular values of each TV for these three decoherence channels, and 
substituting them into Eq. ([9]), we finally obtain the GMQD as follows 



I^Idc = \ [s' {cl + c^) + / + (P' + c,s'? - max {(scO^ , {sc,f ^ + (p^ + c,s'f]] 



4 L 

-^DPC ~ ^ 



s'^{c\ + C2) + C3 — max 
s'^fc? + Co + Co) — max 



[s ci) , [s C2) ,c; 

(Ao^(.%)^(.%)^} 



(23) 
(24) 
(25) 



45| 



For comparison with the dynamics of entanglement, we use the concurrence defined as 4J, 



C = max(0, Ai — A2 — A3 — A4), 



(26) 



where Aj are the square roots of the eigenvalues in descending order of the matrix product 
Q = p {ay ® ay) p* {ay (g) ay) with p* the complex conjugate of the two-qubit density matrix 
P- 

Now we assume p{t) is a smooth function about time t, so we investigate the evolution 
of the correlation quantities (the GMQD, quantum discord and concurrence) along with 
p{t) instead of t. To investigate the behaviors of the GMQD, the quantum discord and the 
concurrence, we consider some examples. The first one is the initial state (ci = 1, C2 = 



C3, C3 = 0.6) [l2| under the PDC. For this case, we obtain 



D^DC(P) = min{Df(p),Df(p)}, 
Df{p) = ^(1 -P)'' Dl{p) = [1 + (1 • (27) 




FIG. 1: (Color online) Discord D{p), concurrence C{p) and the geometry measure of discord D^{p) 
under PDC, as functions of p, for the input state taken as ci = 1, C2 = — C3 and C3 = 0.6. 

In Fig. [H we show that the GMQD is monotonically non-decreasing with respect to the 
quantum discord. When p < 1 — a/3/5, the GMQD decreases while the quantum discord 
keeps constant. It is remarkable that regime where the quantum discord is unaffected by the 
noisy environment is important for the implementation of a quantum computer 12|. The 
entanglement ma y d isappear completely after a finite time, known as entanglement-sudden- 
death (ESD) 30|-l32j. In cases where ESD occurs, contrarily, the quantum discord is more 
robust than the concurrence jisj, so does the GMQD, see Fig. [H The discontinuity of the 
decay rates occurs at p = 1 — a/3/5 (see the dotted line in Fig. [1]). For the GMQD, this 
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kind of sudden change occurs when the maximum singular value of the matrix TV jumps 
from one family to another one, see Df[p) and (p) in Fig. [H 




0.15 




0.05 



FIG. 2: (Color online) Discord D{p), concurrence C{p) and the geometry measure of discord D^{p) 
as functions of p. The left subfigure is plotted for the case of the Bell state (ci = C2 = C3 = — 1) 
under the ADC. The right subfigure is plotted for the case of the states ([29|) with ci = 0.5, 
C2 = 0, C3 = 0.5 and d = —0.5 under the PDC. In the inset of each figure, we plot the values 
of 9 which maximize I{p{p)\{Ei:{9, (p)}) at different p, where Ek{9,(j)) = (1 + n • (t)/2 with n = 
(sin 9 cos (f), sin 9 sin 0, cos 9)^ . 



However, the sudden change in decay rates of the GMQD may not imply that of the 
quantum discord, and vice versa. For instance of the former, we consider the second example, 
a Bell state (ci = 02 = 03 = —1) under the ADC. Substituting q = —1 {i = 1,2,3) into 
Eq. fl2^ . we obtain 



D 



ADC(P) 



mm{Dl{p),Dl{p)} 
Dfip) = 1(1 - 3p + 3p2), Dlip) = 1(1- pf. 



(28) 



From the left subfigure of Fig. 121 we can see that the sudden change in the decay rates of 
the GMQD occurs at p = 0.5, where the quantum discord obtained numerically does not 
display any discontinuity in the first derivative. To demonstrate that the sudden change 
is decay rates of the quantum discord may not imply the GMQD, we consider the third 
example — the input state 

1 



3 




QCTj (g) o"i + doi (g) 1 + (il (g) 0-3 



(29) 



under the PDC Here the parameters are choose as ci = 0.5, C2 = 0, C3 = 0.5 and d = —0.5. 
After some algebras similar to the case where the input states are Bell diagonal states, the 
GMQD is obtained as Dpj^^ = (1 — pY/lG. From the right subfigure of Fig. [21 we can see 
that the sudden change in decay rates of the quantum discord occurs at p ~ 0.22, where the 
GMQD does not display any discontinuity in the first derivative. 

In the following, we give a general analysis on the sudden change in decay rates of the 
GMQD and the quantum discord. Hereafter, we assume that the elements of the two-qubit 
density matrix are smooth functions of p, then the sudden changes in decay rate of both 
the GMQD and the quantum discord are induced by the optimization procedure involved 
in their definitions. For the GMQD, the optimization procedure is to find the closest one x 
among the zero-discord state x = Pi^i ® Pi + P2^2 ® P2- Here Hi and H2 are orthogonal 
projective operators and can be expressed by Hi = (1 + Yl'i=i CjCrj)/2 and H2 = 1 — Hi with 
Cj the components of a unit vector e = (ei, 62, 63)"^ on the Bloch sphere. In Ref. [l9|, Dakic 
et al. obtain the results ^ where fcmax can be expressed by 

fcmax (p) = max [e^K{p) e] . (30) 

|e|=l 

So the sharp features of the GMQD are caused by the optimization over e. Then the closest 
zero-discord state is given by x = pitii pi +^262 ® P2 where tii{p) = (1 + Y^i=i ei(p)o"i)/2 
with e{p) the eigenvector of K{p) with the largest eigenvalue. We do not care about the 
other parameters in the x since they have nothing to do with the GMQD. Hence, the 
sudden change in decay rates of the GMQD corresponds to the sudden change of e{p). 
On the other hand, for the quantum discord, the optimization procedure is to find an 
optimal projective measurement {Ek} to access the classical correlation over the set of 
projective measurement {Ek}. Then the sudden change in the decay rates of the quantum 
discord is induced by the sudden change of {Ek{p)} with respect to p. In a similar way to 
Hj, the optimal projective operator can be represented by -Ei(p) = (1 + X]i=i and 
i?2 = 1 — -El with Ui the components of a unit vector h = (ni, h2, ns)^ on the Bloch sphere, 
or equivalently n = (sin 6' cos 0, sin ^ sin 0, cos 6')-^. In the inset of Fig. |2l we plot 6 for the 
optimized measurement Ei{6,(f)) for different point p. For the second example (a Bell state 
under the ADC), the output state are invariant under a rotation Rz{^p) = exp{iipaf /2) ® 
exp{iipaf /2) and X{p\{Ek}) are invariant under local unitary transformation, so if {Ek{9, 0)} 
is the optimal measurement, then {exp{i{pa^ /2) Ek{9 , 0) exp{—iipa^/2)} = {E{9, (p — ip)} is 
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also the optimal measurement. In other words, the optimization procedure is only relevant 
to 9. In the inset of the left part of Fig. |2l we can see that there is no sudden change of 9. For 
the third example, in the inset of the right part of Fig. [2], we can see that a sudden change 
of 9 occurs at p ~ 0.22, where the decay rate of the quantum discord is discontinuous. And 
we do not care about because the sudden change of the optimal measurement has already 
been reflected through the discontinuity of 9. 

Notice that the optimal Ek{p) and tlk{p) are both projective operators. For general states, 
{Ek{p)} is not the same as {Ilk{p)}, which is reflected in the disagreement of their individual 
sudden change in decay rate. However, at least for the Bell diagonal states, {Ek{p)} is the 
same as itlk{p)}, or equivalently n{p) = e{p). This can be seen from the following analysis. 
In Ref. a], Luo solved the optimization analytically for the Bell diagonal states, and n is 
found to be the eigenvector of RR^ with the largest eigenvalue. On the other hand, in 



Ref. 19| , Dakic et al. showed that e is the eigenvector oi K = xx^ + RR^ with the largest 
eigenvalue. For the Bell states K = RR^ due to x = 0. Hence, we get h{j)) = e{p). So it 
is concluded the sudden change in decay rates of the GMQD that of the quantum discord 
occurs simultaneously if the states are the Bell diagonal states. 



IV. CONCLUSION AND DISCUSSION 



In conclusion, we have considered the dynamics of the GMQD under decoherence. We 
showed that the GMQD of a two-qubit state can be obtained through the singular values 
of a special 3x4 matrix whose elements are the expectation values of the Pauli matrices of 
the two qubits. With the help of the Heisenberg picture, we got the analytic results of the 
geometric measure of the quantum discord for states under three typical kinds of quantum 
decoherence channels. We showed that the sudden change in decay rates of the GMQD does 
not always imply that of the quantum discord, and vice versa. And at least for the Bell 
diagonal states, their individual sudden changes in decay rate are accordance. 

In the present work, we adopt the Hilbert-Schmidt norm as the distance between two 
states. Besides, there exist other quantities for measuring the distance, e.g. the relative 
entropy [20] and the Bures distance. For the Hilbert-Schmidt distance, we show the dis- 
agreement between the GMQD and the quantum discord on reflecting the sudden changes 
in decay rates, so what about other kinds of the distance? Is this disagreement a property of 
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the general geometric measure of quantum discord, or just of the geometry measure based 
on the Hilbert-Schmidt distance? Most recently, the phenomena of the sudden change in 
decay rates have been tested experimentally [18 1. 
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